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ABSTRACT 

We present results from an 850-@m survey of the ~ 5 Myr old A Orionis star-forming region. We 
used the SCUBA-2 camera on the James Clerk Maxwell Telescope to survey a ~ 0?5-diameter circular 
region containing 36 (out of 59) cluster members with infrared excesses indicative of circumstellar disks. 
We detected only one object at > 3 a significance, the Herbig Ae star HD 245185, with a flux density 
of ~ 74 mJy beam -1 corresponding to a dust mass of ~ 150 M®. Stacking the individually undetected 
sources did not produce a significant mean signal but gives an upper limit on the average dust mass 
for A Orionis disks of ~ 3 M$. Our follow-up observations of HD 245185 with the Submillimeter 
Array found weak CO 2-1 line emission with an integrated flux of ~ 170 mJy km s -1 but no 13 CO 
or C 18 0 isotopologue emission at 30 mJy km s -1 sensitivity, suggesting a gas mass of < 1 Mj up . 
The implied gas-to-dust ratio is thus > 50 times lower than the canonical interstellar medium value, 
setting HD 245185 apart from other Herbig Ae disks of similar age, which have been found to be 
gas rich; as HD 245185 also shows signs of accretion, we may be catching it in the final phases of 
disk clearing. Our study of the A Orionis cluster places quantitative constraints on planet formation 
timescales, indicating that at ~ 5 Myr the average disk no longer has sufficient dust and gas to form 
giant planets and perhaps even super Earths; the bulk material has been mostly dispersed or is locked 
in pebbles/planetesimals larger than a few mm in size. 


1. INTRODUCTION 

Circumstellar disks are ubiquitous around proto¬ 
stars, being natural consequences of collapsing molecular 
clouds conserving angular momentum while contracting 
over two orders of magnitude in size. Space-based in¬ 
frared (IR) surveys of young stellar clusters have placed 
strong constraints on disk lifetimes, revealing an expo¬ 
nentially decreasing disk occurrence with age and a char¬ 
acteristic disk lifetime of only a few Myr (e.g., Ribas et al. 
2014; Mamajek 2009; Hernandez et al. 2007). However, 
IR observations are sensitive to small amounts of dust 
due to optical depth effects, saturating at less than a 
lunar mass of material spread over ~ 100 AU, making 
IR excess a poor tracer of total dust mass and thus the 
planet-forming capacity of disks. 

Instead, sub-mm emission is used to probe total disk 
mass. Continuum emission at these longer wavelengths is 
optically thin and therefore can be directly related to the 
total mass of the emitting dust (Hildebrand 1983). Sur¬ 
veys of young stellar clusters are inconclusive on whether 
sub-mm emission systematically declines with age, which 
would reflect steady disk dispersal and grain growth. An¬ 
drews et al. (2013) showed that, when accounting for 
different stellar populations and survey sensitivities, the 
disk populations in Ophiuclius (~ 1 Myr; Andrews & 
Williams 2007) and IC348 (~ 2 Myr; Lee et al. 2011) are 
consistent with that of Taurus (~ 2 Myr; Andrews et al. 

2013) , while the discrepancy with the evolved Upper Sco 
region (~ 10 Myr; Mathews et al. 2012; Carpenter et al. 

2014) is only marginal (2.5 ct). However, Williams et al. 
(2013) has also shown that the disk population in the 
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intermediate-aged a Orionis region (~ 4 Myr) is distin¬ 
guishable from that of Taurus at > 3cr significance. 

It is important to note that estimating total disk mass 
(gas + dust) from sub-mm continuum flux, which probes 
only the dust, requires assuming a gas-to-dust ratio. The 
canonical interstellar medium (ISM) gas-to-dust ratio of 
~ 100 (Bohlin et al. 1978) is often applied, such that 
Afdisfc ~ 100 Mdust ■ This has been a standard approach 
as obtaining accurate gas masses is hindered by weak 
line emission as well as complicated disk chemistry and 
radiative transfer. Although Class 0/1 disks around em¬ 
bedded protostars presumably form with an inherited 
ISM gas-to-dust ratio, they evolve rapidly to Class II 
protoplanetary disks around optically visible stars, and 
eventually to Class III debris disks with negligible gas. 
Thus this extrapolation of two orders of magnitude is 
not only large but also uncertain. Using a recently de¬ 
veloped technique to independently measure gas masses 
from CO isotopologues, Williams & Best (2014) found 
that the gas-to-dust ratios of nine protoplanetary disks 
in the young (~ 2 Myr) Taurus association are already 
an order of magnitude lower than the ISM value and may 
vary significantly from disk to disk. This suggests that 
assuming a constant ISM gas-to-dust ratio, for even very 
young disks, results in significantly over-estimated total 
disk masses. 

In this paper we focus on the A Orionis cluster in or¬ 
der to further understanding of disk mass evolution and 
better constrain planet formation timescales. A Orionis 
is an evolved (~ 5 Myr; Dolan & Mathieu 2001) star¬ 
forming region, centered on a tight core of OB stars 
and surrounded by a 30-pc radius ring of dense dust 
and gas, possibly formed by a supernova that occurred 
~ 1 Myr ago (e.g., Dolan & Mathieu 1999). For our 
survey we used the Sub-mm Common User Bolome¬ 
ter Array (SCUBA-2) camera on the 15-m James Clerk 
Maxwell Telescope (JCMT) atop Maunakea, as it com- 
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Table 1 

Observing log 


Date hnt a T 2 25 b 


JCMT 


2013/09/23 

1.5 

0.10 

2014/01/17 

2.2 

0.20 

2014/01/18 

1.7 

0.17 

2014/01/20 

0.4 

0.16 

2014/08/23 

1.3 

0.09 

2014/08/26 

1.3 

0.09 

2014/10/13 

2.2 

0.09 

2014/11/12 

1.8 

0.13 

2015/01/11 

3.9 

0.11 

2015/04/17 

1.4 

0.05 

2015/04/19 

1.8 

0.04 


SMA 


2015/02/02 4.6 0.14 

a On-source integration time in hrs. 

k Average zenith optical depth at 225 GHz. 

bines high sensitivity with large fields of view, enabling 
unbiased surveys of young stellar clusters on compara¬ 
ble scales to the aforementioned space-based IR surveys. 
Our study focuses on the area < 0.5° from the central 
core, which encompasses the several hundred pre-main 
sequence (PMS) stars typically associated with the A 
Orionis cluster. Although the global Initial Mass Func¬ 
tion (IMF) of A Orionis is similar to the field, this inner 
region is marginally overpopulated by OB stars (Dolan 
& Mathieu 2001). The region is too distant for reliable 
Hipparcos parallax measurements, thus we adopt the dis¬ 
tance of ~ 450 pc derived from Stromgren photometry 
of the OB stars (Dolan & Mathieu 2001). Reddening 
toward A Orionis is low at E(B — V) ~ 0.12 (Diplas & 
Savage 1994) and thus extinction effects should be neg¬ 
ligible. See Mathieu (2008) for a detailed account of the 
cluster properties. 

We begin by describing our sample selection and ob¬ 
servations in Section 2, then present the survey results in 
Section 3. In Section 4, we compare our results to other 
young clusters and take a detailed look at our only de¬ 
tection (HD 245185) in order to analyze the implications 
for disk evolution and planet formation. In Section 5, we 
summarize our paper and discuss areas of future work. 

2. SAMPLE & OBSERVATIONS 
2 . 1 . Sample selection 

We formed our sample of disk-bearing A Orionis mem¬ 
bers by combining results from Hernandez et al. (2009) 
and Hernandez et al. (2010). They surveyed the region 
with Spitzer to identify cluster members and classify their 
circumstellar disks according to IR excess. Hernandez 
et al. (2009) studied the intermediate-mass population, 
finding 29 members earlier than F5 but only 10 bearing 
disks. Hernandez et al. (2010) studied the lower-mass 
population, finding 436 members down to the substel- 
lar limit but only 49 with disks. This low disk fraction 
is consistent with the ~ 5 Myr age of the cluster (e.g., 
Ribas et al. 2014; Mamajek 2009; Hernandez et al. 2007). 

Hernandez et al. (2009) classified the disks around 
intermediate-mass members according to their excess 
emission in Spitzer IRAC/MIPS bands. They classified 



Figure 1. The histogram of flux densities for all pixels within 
the usable region of the SCUBA-2 map (Figure 2). The one 
detection (HD 245185) was masked out as applying a matched 
filter (Section 2.2.1) creates negative troughs around strong 
sources, resulting in a non-Gaussian tail of negative flux den¬ 
sities. The Gaussian fit to the histogram (black line) has 
a dispersion of 2.9 mjy beam -1 . The surplus at high flux 
densities shows the small number of possible detections. 

nine sources with moderate excess as “Debris” disks and 
one source with large excess in all Spitzer bands as an op¬ 
tically thick “Primordial” disk. Hernandez et al. (2010) 
classified disks around lower-mass members according to 
their spectral energy distribution (SED) IR slopes. They 
classified optically thick disks with the largest excesses as 
“Thk,” evolved disks with smaller excesses as “Ev,” and 
(pre-) transition disks with signs of inner disk clearing as 
“PTD” or “TD.” For simplicity, we group these classifi¬ 
cations into optically thick disks (Primordial, Thk) and 
optically thin disks (TD, PTD, Ev, Debris). 

2 . 2 . Observations 
2.2.1. SCUBA-2 

We observed a ~ 0? 5-diameter circular region toward 
the A Orionis cluster with SCUBA-2 for a total of ~ 19.2 
hrs. The observations, summarized in Table 1, were 
taken in queue mode over nine nights from Septem¬ 
ber 2013 to April 2015 (program IDs: M13BH12A, 
M14AH16A, M14BH050, M15AH01). JCMT defines 
weather condition bands based on zenith optical depths 
at 225 GHz (t 2 25 ); our data were taken in JCMT band 
2 (0.05 < r 22 5 < 0.08), band 3 (0.08 < r 225 < 0.12), or 
band 4 (0.12 < T 225 < 0.20) weather. We chose the point¬ 
ing center of a = 5 h 35 m 4 s and 5 = +9°56 , 60" to maxi¬ 
mize the number of our targets within the map area. The 
PONG1800 mapping mode was used to provide nearly uni¬ 
form sensitivity over the large circular field. Flats were 
taken before and after each scan. 

We reduced the data using the Dynamic Iterative Map 
Maker (DIMM) in the STARLINK/smurf software pack¬ 
age (Chapin et al. 2013). The DIMM performs pre¬ 
processing and cleaning of the raw data (e.g., dark sub¬ 
traction, concatenation, flat fielding) to create reduced 
science maps with a default pixel scale of 4" at 850 
fim. These reduced maps are then flux calibrated and 
co-added using the Pipeline for Combining and Analyz¬ 
ing Reduced Data (PICARD), a post-processing software 
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Figure 2. Our SCUBA-2 850-/im map of A Orionis, showing locations of disk-bearing members covered by the usable map area 
(Section 2.2.1). Squares indicate optically thick disks and circles indicate optically thin disks (Section 2.1). The only significant 
detection, HD 245185, is highlighted in red. The blue diamonds indicate possible sub-mm galaxies (Section 3.3 & Table 4). 


also included in the STARLINK/SMURF package. 

We reduced each scan individually. We first ran DIMM 
on the individual chips for a given scan (SCUBA-2 is a 
mosaic of four detector chips per waveband). We ap¬ 
plied a “blank field” DIMM configuration file, which is 
commonly used for deep-field extragalactic surveys to en¬ 
hance faint point source detection (circumstellar disks 
at the distance of A Orionis should be point sources at 
the 15" resolution of SCUBA-2 at 850 //m). We then 
used PICARD to co-add the reduced chip maps into a 
single scan map. To maximize point-source detectabil¬ 
ity, we also applied a matched filter, which suppresses 
residual large-scale noise by self-subtracting a Gaussian- 
smoothed version of the map. 

We then flux calibrated each scan. The flux density 
scale was determined from regular observations of bright 
point sources, either Uranus or the protoplanetary neb¬ 
ula CRL 618. The calibrators were reduced using the 
same method described above for the science targets in 
order to derive a flux calibration factor (FCF), which 
scales data from units of pW to Jy. Several different 
FCFs can be applied depending on the data and science 
needs. Because we applied a matched filter to our science 


maps, we used FCF matc h, which sets the peak pixel to the 
total flux density of the point source. FCF match is very 
stable against beam distortion; we found FCF matc h = 
672 ± 36 Jy pW -1 beam -1 across our observations, and 
this ~ 5% variation represents the absolute flux calibra¬ 
tion uncertainty. We applied an FCF matc h to each re¬ 
duced scan map, using the calibrator observed closest in 
time. Before applying the calibration, we inspected the 
integrity of each calibrator by re-reducing them using the 
“bright compact” DIMM configuration file and checking 
that the resulting FCF pea k values matched the nominal 
values from Dempsey et al. (2013) to within 5%. 

After each scan was reduced and flux calibrated, we 
used PICARD to combine all the scans for a given night, 
then combined all the nightly maps into our final 850- 
(im map, shown in Figure 2. The usable map area was 
determined by iteratively cropping the noisy outer map 
regions until the histogram of pixel flux densities repre¬ 
sented a Gaussian distribution (Figure 1). This resulted 
in a ~ 0?45-diameter map with a median rms noise of 
~ 2.9 mJy beam -1 . The map contains 36 of the 59 
known A Orionis members with circumstellar disks (Sec¬ 
tion 2.1). These sources are listed in Table 2. 
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Although SCUBA-2 simultaneously images both 850- 
and 450-/im wavebands, the weather conditions were too 
poor for a 450-/mr survey. However, we recovered our 
one 850-/xm detection (HD 245185; Section 4.3) at this 
shorter wavelength and thus reduced the 450-/im map in 
a similar manner as for the 850-/im map described above. 

2.2.2. Submillimeter Array 

We followed up our only SCUBA-2 detection (HD 
245185; Section 4.3) with interferometric observations 
using the Submillimeter Array (SMA). We used the SMA 
1300-/im band with the correlator configured to also 
cover CO 2 1 line emission at 230.538 GHz. This al¬ 
lowed us to measure the continuum flux at a slightly 
longer wavelength to extend the SED (Figure 5) as well 
as search for molecular gas in the disk. 

Data were taken on 2015-02-02 in extended configu¬ 
ration with baselines of ~ 30-200 m. Weather condi¬ 
tions were average with precipitable water vapor of ~ 2 
mm. The observational setup was typical for SMA, con¬ 
sisting of a long integration on a passband calibrator 
(3C279) followed by 20-min integrations on the source 
interleaved with 3-5-min integrations on two gain cali¬ 
brators (0510+180 and 3C120). The absolute flux scale 
was derived from observations of Ganymede and is ac¬ 
curate to ~ 20%. The data reduction followed standard 
procedures, using the MIR software package to calibrate 
the visibilities and MIRIAD for final imaging. The beam 
sizes of the final maps were l'.'OO x 0'/86 for the 1300-/mi 
continuum and 0 , /90 x 0'/82 for the CO 2-1 line emission. 

3. RESULTS 
3.1. SCUBA-2 

For each source, we measured the total flux density 
using the pixel value at its Two Micron All-Sky Survey 
(2MASS; Skrutskie et al. 2006) source location. We took 
this approach because we calibrated our science maps us¬ 
ing FCF matc h (Section 2.2.1), which sets the peak pixel 
value to the total flux density for point sources. We 
took the value at the 2MASS source location, rather 
than searching for the peak pixel within the beam, as 
the SCUBA-2 850-/im pixel size (4") is larger than the 
combined 2MASS positional errors (< 1") and JCMT 
nominal pointing performance (1.5-2.0"). We estimated 
ler statistical uncertainties from the rms noise within an 
annulus between 20" and 60" in radius. 

The results are tabulated in Table 2. Only one source, 
the Herbig Ae star HD 245185, was detected with > 3 ct 
significance (F 850 = 74.07 + 2.63 mJy). We were also able 
to detected HD 245185 at 450 /mi to ~ 6a significance 
(F 450 = 326 ± 52 mJy). These are the first sub-mm ob¬ 
servations of HD 245185, although the target is a known 
mm source (Section 3.2). 

We stacked the individually non-detected sources on 
their 2MASS positions, then measured the stacked flux 
(Fstack) and associated statistical error (astack) in the 
same manner as described above for individual sources in 
order to check for a statistically significant mean signal. 
As reported in Table 3, we did not find a significant mean 
signal in the stacked image, even when considering only 
non-detections with optically thick disks. We verified 
this by calculating the mean flux density ( F ),) and stan¬ 
dard error on the mean (cr^) for all the non-detections, as 


Table 2 

SCUBA-2 850-/im Survey Results 


ID a 

RA b 

Dec b 

SpT c 

F 85 o d 

e^so 11 

Class a 

1624 

83.5207 

9.9511 

M7.0 

4.57 

3.39 

Thk 

2088 

83.5803 

9.8076 

M7.0 

3.92 

2.71 

Thk 

2357 

83.6163 

10.0467 


-2.44 

2.93 

Ev 

2404 

83.6215 

9.7853 


3.00 

2.83 

Ev 

2712 

83.6597 

9.9699 

M5.0 

2.70 

2.50 

Ev 

2989 

83.6926 

9.9271 

M5.5 

-2.99 

2.40 

Thk 

2993 

83.6930 

10.0422 

M5.0 

-0.11 

2.88 

Ev 

3360 

83.7350 

9.9179 

M4.0 

4.77 

2.44 

Thk 

3506 

83.7507 

9.8780 


-0.96 

2.83 

Ev 

3597 

83.7614 

9.9466 

M4.0 

-2.09 

2.42 

Ev 

3710 

83.7762 

9.7816 

M5.5 

-1.80 

2.91 

Ev 

3746 

83.7795 

9.9004 


-0.62 

2.55 

Thk 

3785 

83.7847 

9.7149 


-6.35 

3.43 

Ev 

3887 

83.7964 

9.9554 

M5.0 

-1.87 

2.95 

Thk 

3942 

83.8023 

9.8864 

M3.0 

-2.95 

2.35 

Ev 

4021 

83.8139 

9.8103 

M3.0 

2.90 

2.73 

PTD 

4111 

83.8247 

9.9492 


3.35 

2.52 

TD 

4126 

83.8258 

9.8734 

M2.0 

-2.89 

2.13 

Thk 

4155 

83.8294 

9.9153 

K2 

6.11 

2.63 

Ev 

4163 

83.8298 

9.9118 

M4.0 

2.17 

2.54 

Thk 

4187 

83.8330 

10.0435 

M1.0 

5.17 

3.43 

Thk 

4255 

83.8397 

9.8914 


-3.74 

2.31 

Thk 

4363 

83.8517 

9.8977 

M5.5 

-1.25 

2.79 

Thk 

4407 

83.8557 

10.1440 


-1.73 

2.95 

Thk 

4520 

83.8686 

10.0410 

M3.5 

2.48 

2.83 

Thk 

4531 

83.8699 

9.9028 

M5.5 

1.11 

2.64 

Thk 

4817 

83.9051 

9.9477 


3.55 

2.89 

Thk 

4916 

83.9160 

9.8900 

M6.5 

-0.67 

2.72 

Thk 

5042 

83.9348 

10.0984 

M6.0 

-0.41 

2.96 

Thk 

5447 

83.9827 

9.9394 

M3.5 

-3.56 

3.09 

Thk 

7951 

83.7831 

10.0017 

M6.0 

-1.14 

2.58 

Thk 

7957 

83.6512 

9.9256 

M8.0 

2.31 

2.72 

Thk 

HD36894 

83.8196 

9.7776 

B2 

-0.58 

2.90 

Debris 

HD245168 

83.7624 

9.9345 

B9 

-1.36 

2.47 

Debris 

HD245185 

83.7900 

10.0310 

A0 

74.07 

2.63 

Primordial 

HD245275 

83.9369 

9.9234 

A5 

5.27 

2.72 

Debris 


a From Hernandez et al. (2009) and Hernandez et al. (2010). 
b Units of deg (J2000). 

c Spectral types from literature (Sacco et al. 2008; Barrado y 
Navascues et al. 2007; Bayo et al. 2011; Franciosini & Sacco 2011) 
d Units of mJy beam -1 . 


Table 3 

Stacking analysis 


Type 

N a 

rp b 

r stack 

&stack 

cr 

°> b 

Thick 

19 

0.46 

0.72 

0.18 

0.60 

All 

35 

0.40 

0.46 

0.30 

0.46 


a Number of non-detected sources used. 
° Units of mJy beam -1 . 


well as just the non-detections with optically thick disks, 
the results of which are also listed in Table 3. 

3.2. Submillimeter Array 

Figure 3 shows the results from our SMA observations. 
We clearly detected 1300-/im continuum emission from 
HD 245185 where +1300 =30 + 6 mJy. HD 245185 has 
been previously detected at mm wavelengths: Mannings 
& Sargent (1997) reported F 1300 = 44 + 12 mJy and 
+2600 = 6.5 ± 1.2, which is consistent with our results. 
The source was unresolved in the l'/OO x 0"86 SMA beam, 
thereby constraining the disk radius to < 200 AU, which 
is an improvement on the previous upper limit of ~ 680 
AU from Mannings & Sargent (1997). 
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Table 4 

Candidate sub-mm galaxies 


Name 

RA a 

Dec a 

Fsso h 

&850 b 

SMG J053450.5-100120 

83.7106 

10.0222 

13.23 

2.63 

SMG J053552.0-100304 

83.9667 

10.0510 

13.57 

2.84 

SMG J053550.4-100443 

83.9600 

10.0787 

13.30 

2.61 


a Units of deg (J2000). 

° Units of mjy beam -1 . 


used a 4.5cr threshold for detection of possible SMGs. We 
found 3 candidate SMGs that lacked counterparts within 
15" at all other wavelengths. This matches expectations 
from the number counts model of Chen et al. (2013b) 
and is the same number found by Williams et al. (2013) 
in their SCUBA-2 map of er Orionis with similar size and 
sensitivity. The locations and fluxes of these candidate 
SMGs are given in Table 4. 



Figure 3. Results from our SMA observations of HD 245185 
(Section 3.2). The top panel shows the strong but unresolved 
1300-/rm continuum emission where F 1300 = 30 ± 6 mjy 
(contours in units of mjy). The bottom panel shows the CO 
2- 1 line emission at 230.538 GHz with a weak integrated 
flux of 170 ± 30 mjy km s -1 . A fitted Gaussian (dashed 
red curve) is centered at 11.0 ± 0.1 km s -1 with a narrow 
velocity profile of FWHM = 0.3 ± 0.1 km s -1 . 


Weak CO 2-1 line emission was found at u ra( j = 11.0 ± 
0.1 km s -1 with a peak of 0.31 ± 0.05 Jy and a velocity 
resolution of 0.1 km s -1 . The velocity profile is very 
narrow with FWHM = 0.3 ± 0.1 km s -1 , suggesting a 
face-on disk geometry. The integrated line flux is weak 
at 170 ±30 mjy km s _1 . No 13 CO or C 18 0 emission was 
detected at the rms noise level of 30 mjy km s _1 . 

3.3. Background sources 

Sensitive, large-scale SCUBA-2 surveys similar to this 
work are used to search for faint sub-mm galaxies 
(SMGs). These dusty, star-forming galaxies at high red- 
shift are interesting because they are thought to have 
been instrumental in the build-up of massive galaxies 
(Casey et al. 2014). Based on our SCUBA-2 map size, we 
expected one 4cr noise spike (Chen et al. 2013a) and thus 


4. DISCUSSION 
4.1. Dust Masses 

We can estimate the dust mass of each disk using the 
850-/im continuum fluxes in Table 2. Assuming the dust 
emission at sub-mm wavelengths is optically thin, the 
sub-mm continuum flux (F v ) is directly related to the 
dust mass (M dust ) as in Hildebrand (1983): 


Afdust 


F„d 2 

K V B u (r Anst y 


( 1 ) 


where B v is the Planck function and we assume a char¬ 
acteristic dust temperature of T Aust = 20 K, which is 
the median found for Taurus disks (Andrews & Williams 
2005). The dust grain opacity, k„, is taken as 10 cm 2 
g -1 at 1000 GHz with an opacity power-law index /3 = 1 
(Beckwith et al. 1990). Using the distance to A Orionis 
of d « 450 pc (Dolan & Mathieu 2001), we can estimate 
dust masses from F 850 (mjy) using: 


Miust ~ 5-8 x 10 6 F850 M 0 . (2) 

This gives a dust mass for the disk around HD 245185 
of ~ 150 M ffi and a 3er upper limit on the average dust 
mass for all other A Orionis disks of ~ 3 Mg. 


4.2. Comparison with other star-forming regions 

In Figure 4, we compare the dust mass distribution of 
optically thick Class II disks in A Orionis with those of 
several other star-forming regions taken from the litera¬ 
ture with ages spanning ~ 1-10 Myr. We compare only 
Class II disks as they can still exhibit relatively strong 
sub-mm emission but lack confusion from large circum- 
stcllar envelopes dominating Class 0/1 disks. We com¬ 
pare only dust masses, rather than total disk masses, as 
gas-to-dust ratios are highly uncertain and may vary sig¬ 
nificantly from disk to disk, even within the same cluster 
(e.g., Williams & Best 2014). As slightly different meth¬ 
ods are often applied to derive dust masses, we did not 
take the values directly from the literature. Rather, we 
translated the published sub-mm continuum fluxes into 
dust masses using versions of Equation 1 scaled to the 
distances of the clusters and the observation wavelengths 
of the surveys. The survey sensitivities (shown by arrows 
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log (Mdust) (M®) 

Figure 4. Dust mass distributions for Class II disks in clus¬ 
ters with ages spanning ~ 1-10 Myr: Taurus (Andrews & 
Williams 2005), cr Orionis (Williams et al. 2013), A Orionis 
(this work), and Upper Sco (Carpenter et al. 2014). The ar¬ 
rows indicate 3<r survey sensitivities and the cluster ages (in 
Myr) are shown next to the cluster name. 

in Figure 4) are determined from their respective median 
3 a statistical errors. 

The surveys in Figure 4 are all sensitive to dust masses 
> 30 M®. This allows us to study how the planet-forming 
capacity of disks evolves with time, assuming gas-to-dust 
ratios of ~ 2-40 (e.g., Williams & Best 2014). From 
Figure 4, it appears that over ~ 1- 10 Myr timescales 
there is substantial evolution in the frequency of disks 
capable of forming systems with multiple giant planets: 
such massive disks drop in frequency by roughly an or¬ 
der of magnitude by ~ 5 Myr, presumably because the 
disk material has mostly dispersed or is locked in peb- 
bles/planetesimals > 1 mm in size. During ~ 5-10 Myr, 
only outlier disks maintain such high masses. Note that 
a proper statistical comparison between disk populations 
should account for different stellar mass distributions, in¬ 
dividual disk sensitivities, and survey incompletenesses 
(e.g., Andrews et al. 2013); however, this approach was 
impractical here due to the single detection in A Orionis. 

Figure 4 also highlights the particularly high dust mass 
of HD 245185 for its age, making it an outlier in the 
A Orionis cluster—and perhaps the disk population in 
general. Indeed, HD 245185 resides far from our survey 
sensitivity limit and we did not detect any other disks 
in this region. Moreover, Hernandez et al. (2009) noted 
that the 24-/im excess from HD 245185 was over twice as 
large as other disks around intermediate-mass stars with 
similar ages. Ribas et al. (2015) also showed that proto¬ 
planetary disks represent only 2% of the disk population 
at ~ 5 Myr for stars > 2 Mg. Below we discuss possible 
physical reasons for this outlier status, namely cluster 
non-membership and gap-opening planet formation. 

4.3. HD 245185 

The only significant detection in our 850-/xm survey of 
A Orionis was HD 245185, a well-known Herbig Ae star 
of ~ 2.5 Mg (e.g., Folsom et al. 2012). With a detec¬ 
tion significance of ~ 28cr, its dust mass is > 50 times 



logA (pm) 

Figure 5. The optical-mm wavelength SED for HD 245185. 
Errors are all smaller than the symbol sizes. The gray line is 
a NextGen model photosphere (Hauschildt et al. 1999) for a 
9600 K star with log g = 4.0 and [M/H] = 0 (Folsom et al. 
2012), scaled to a distance of 450 pc. Extinction corrections 
are not included due to the low reddening toward A Orionis. 

larger than the average A Orionis disk. Its strong sub- 
nmi continuum emission, in addition to its large Spitzer 
IR excesses (Hernandez et al. 2009), suggests that HD 
245185 hosts a massive primordial disk. However, proto¬ 
planetary disks around intermediate-mass stars are sup¬ 
posed to disperse twice as fast compared to low-mass 
stars (Ribas et al. 2015), raising the question of why 
such a massive primordial disk has persisted around HD 
245185. In this section we present a more detailed anal¬ 
ysis of HD 245185 to investigate this question. 

4.3.1. Spectral energy distribution 

Figure 5 shows the optical-mm wavelength SED of HD 
245185 (values given in Table 5). The circumstellar disk 
around the central star produces the IR excesses above 
the expected stellar photosphere emission (gray line in 
Figure 5). Hernandez et al. (2009) classified the disk as 
primordial due to its large mid-IR excesses in all Spitzer 
bands. Our SED shows that this excess exists even at 
near-IR wavelengths, which suggests that the disk ex¬ 
tends down to very small radii (< 1 AU). 

However, the mid-IR dip in the SED is a strong indi¬ 
cator of a pre-transition disk, which has evolved past the 
primordial disk phase by developing a gap in the disk 
structure (Espaillat et al. 2014). The large disk mass 
(~ 150 M®; Section 4.1) rules out a highly evolved de¬ 
bris disk as these typically have much lower dust masses 
of ~ 10' 1 M® (see Figure 3 in Wyatt 2008). We also fit 
the tail end of the SED to calculate the sub-mm spectral 
index, assuming F v oc v a at these longer wavelengths. 
We found a shallow slope of a = 2.23 ± 0.04, close to 
the blackbody value suggesting grain growth to at least 
mm-size objects (Draine 2006). This a value is similar 
to those found for other disks (e.g., Testi et al. 2014). 

Note that we show only one WISE band in Figure 5 
for clarity, as the others overlap those of Spitzer. In 
Table 5 we include all four WISE measurements, which 
are consistent with the shape of the SED traced by the 
Spitzer measurements in Figure 5. Thus there is no ev- 
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Table 5 

Spectral energy distribution 


^ef f 

(Axio 3 ) 

F 

(Jy) 

Source [Filter] 


3.660 

0.1725 ± 0.0048 

de Winter et al. 2001 

U] 

4.380 

0.3989 ±0.0110 

de Winter et al. 2001 

B} 

5.450 

0.3878 ±0.0107 

de Winter et al. 2001 

V ] 

6.410 

0.3422 ± 0.0094 

de Winter et al. 2001 

Re] 

7.980 

0.2878 ± 0.0080 

de Winter et al. 2001 

UJ 

12.35 

0.3063 ± 0.0085 

2MASS 



16.62 

0.3197 ± 0.0147 

2MASS 

H] 


21.59 

0.4130 ±0.0152 

2MASS 

A S ] 


33.68 

0.5414 ±0.0299 

WISE [1 



35.50 

0.5983 ±0.0165 

Spitzer [IRAC-1] 


44.93 

0.5945 ± 0.0164 

Spitzer [IRAC-2] 


46.18 

0.6073 ±0.0145 

WISE [2] 


57.31 

0.7269 ±0.0201 

Spitzer [IRAC-3] 


78.72 

1.6945 ± 0.0468 

Spitzer [IRAC-4] 


120.8 

3.2591 ± 0.0360 

WISE [3] 


221.9 

5.2134 ±0.0720 

WISE [4] 


236.8 

4.7021 ± 0.0001 

Spitzer [MIPS 24 /im] 


4500 

0.3262 ± 0.0524 

This work 


8500 

0.0741 ± 0.0042 

This work 


13000 

0.0290 ±0.0058 

This work 


26000 

0.0065 ± 0.0012 

Mannings & Sargent 1997 


idence of HD 245185 exhibiting mid-IR variability over 
a several-yr timescale, which could indicate major colli¬ 
sions occurring in the disk (Meng et al. 2015) 



Figure 6. Archival CFHT/ESPaDOnS spectra of HD 245185. 
Top panel shows [O i] 6300A line emission with a median er¬ 
ror bar given for reference. The rc-axis is centered on the 
stellar radial velocity of v la .d = 11.0 km s -1 (Section 3.2); the 
peak of the Gaussian fit (red line) after masking out telluric 
absorption features (gray regions) is offset by 5 km s -1 but 
within measurement errors. Bottom panel shows Ha 6563A 
line emission, likely from accretion as it is strong, broadened, 
double-peaked, and axisymmetric. 


4.3.2. CO emission 

We detected weak CO 2-1 line emission from HD 
245185 (Section 3.2), indicating the presence of molec¬ 
ular gas in the disk. Converting integrated CO 2-1 line 
flux into a total gas mass is complicated by the optical 
thickness of the line emission. Instead, our best con¬ 
straints come from comparing our upper limits on the 
CO isotopologues to the grid of models in Williams & 
Best (2014), who used these optically thin lines to pre¬ 
dict total gas masses to within a factor of three. Our 3 ct 
upper limits to the 13 CO and C 18 0 line luminosities for 
HD 245185 (~8x 10 4 Jy km s” 1 pc 2 ) only cover model 
disks with M gas < 1 Mj up (see Figure 6 in Williams & 
Best 2014). This implies a total gas mass for HD 245185 
of < 1 Mj up , which when combined with our measured 
dust mass of ~ 150 M® (Section 4.1) gives a notably low 
gas-to-dust ratio of < 2. Although we cannot rule out a 
small, optically thick disk as an alternative explanation 
for the low CO 2-1 line emission (we only have an upper 
limit on the disk radius of ~ 200 AU; Section 3.2), disks 
around Herbig Ae stars are typically > 100 AU in radius. 

This substantial gas depletion sets HD 245185 apart 
from other disks around Herbig Ae stars of similar age, 
which have been found to remain relatively gas rich. For 
example, the transition disk HD 169142 (~ 6 Myr) has 
M gas ~ 20 Mj up (Panic et al. 2008) and Md us t ~ 20 Mj up 
(Raman et al. 2006), thereby maintaining its initial ISM 
gas-to-dust ratio of ~ 100. The massive disk HD 163296 
4 Myr) has M gas ~ 50 M Jup and M dust ~ 0.3 M Jup 
(Williams & Best 2014), suggesting a gas-to-dust ratio of 
~ 170 (Qi et al. 2011 also estimate a ratio of ~ 150). The 
distinction from HD 245185 is significant, even when con¬ 
sidering that gas-to-dust ratios are typically uncertain by 
a factor of ~ 10. 


4.3.3. [Ol] & Ha emission 

To further investigate the gas component of the disk 
around HD 245185, we analyzed the non-thermal [Ol] 
6300A emission line. This forbidden emission arises 
when UV photons from the host star photodissociate OH 
molecules in the surface layers of the circumstellar disk 
to create excited oxygen atoms, some of which are in the 
upper state of the 6300A line (van Dishoeck & Dalgarno 
1984). [Ol] 6300A line luminosities have been empir¬ 
ically correlated to accretion rates for Herbig Ae stars 
(Mendigutfa et al. 2011), which can be physically ex¬ 
plained if the UV excess from the accretion shock is illu¬ 
minating the disk to excite the oxygen atoms. Although 
[O i] 6300A emission is clearly not a good measure of to¬ 
tal gas mass, it can be used to investigate the presence 
of gas in the inner disk of HD 245185 despite its depleted 
overall gas content (Section 4.3.2). 

Acke et al. (2005) previously analyzed the [O i] 6300A 
emission line from HD 245185 using a high-resolution 
spectrum (R = 45,000) but with limited S/N (see their 
Figure 2). We therefore re-visited this analysis with 
higher resolution (i? = 65,000) and higher S/N (~ 80 per 
resolution element at 6300 A) archival spectra from the 
Echclle Spectropolarimetric Device for the Observation 
of Stars (ESPaDOnS) instrument on the Canada-France- 
Hawaii Telescope (CFHT). We median combined the four 
spectra taken on 2005-02-20, normalized the spectrum to 
the fitted continuum, masked out telluric absorption fea¬ 
tures, and fit a Gaussian to the spectral line. The result¬ 
ing spectrum and fit are shown in Figure 6. We found 
|EW| = 0.13±0.08 A and FWHM = 43±5 km s"\ both 
of which are typical for Herbig Ae/Be stars with pure 
emission line profiles (e.g., group I in Acke et al. 2005). 
















The offset of the emission line peak from our derived stel¬ 
lar radial velocity of u ra( j = 11.0 km s -1 (Section 3.2) is 
within measurement errors. 

These results suggest that HD 245185 has the typical 
[O i] 6300A line emission expected from disks around ac¬ 
creting Herbig Ae stars. This accretion signature is sup¬ 
ported by the Ha 6563A line emission from HD 245185, 
also seen in its archival CFHT/ESPaDOnS spectrum 
(Figure 6). The Ha emission line is strong (|EWhq| ~ 34 
A), broadened, double-peaked, and asymmetric. Indeed, 
Manoj et al. (2006) previously found |EWh«| ~ 25 A and 
suggested this as a signature of accretion due to correla¬ 
tions between Ha emission and near-IR excess. Donehew 
& Brittain (2011) also used veiling of the Balmer discon¬ 
tinuity to derive an accretion rate of 0.63 ± 0.41 x 10 - ' 
M 0 yr -1 , which is typical for HD 245185’s stellar mass 
and age (e.g., Mendigutia et al. 2011; Liu et al. 2011). 
This ongoing accretion is surprising given the depleted 
overall gas content for this disk (Section 4.3.2), and sug¬ 
gests that we may be catching HD 245185 at the final 
stages of disk clearing. 

4.3.4. Cluster membership 

HD 245185 is a clear outlier in the A Orionis star¬ 
forming region in terms of its circumstellar disk prop¬ 
erties. This brings into question whether HD 245185 is a 
true member of the cluster. Dynamically, HD 245185 is 
a debatable member of A Orionis. The proper motion of 
HD 245185 is p a cos(5) « 0.67T1.3 and p s « -1.60±1.3 
mas yr -1 (Roser et al. 2008), which is in very good agree¬ 
ment with the average proper motion of the early-type 
members of A Orionis (see Figure 1 in Hernandez et al. 
2009). However, the precise radial velocity derived from 
our SMA observations of u ra( j = 11 ± 0.1 km s -1 (Sec¬ 
tion 3.2) is significantly different from the tight distri¬ 
bution of radial velocities for A Orionis cluster members 
of F ra d = 24.3 ± 2.3 km s -1 (Dolan & Mathieu 2001). 
Peculiar kinematics of high-mass stars in the A Orionis 
cluster have been previously attributed to a supernova 
that may have occurred ~ 1 Myr ago in the A Orio¬ 
nis stellar system (e.g., Dolan & Mathieu 2001; Cunha 
& Smith 1996), suggesting kinematics alone are an un¬ 
reliable criteria for membership. However, HD 245185 
is unlikely to have been affected by such a supernova, 
at least kinematically, as the projected distance between 
HD 245185 and the A Orionis star is almost a parsec. 

Interestingly, the kinematics of HD 245185 are consis¬ 
tent with those of the nearby, younger Taurus region. 
Taurus has u ra d = 16 ± 6 km s -1 (Bertout & Genova 
2006), which is consistent with the ii ra d of HD 245185 
derived in this work. The proper motion of HD 245185 
is also within 3<r of the average proper motion of Tau¬ 
rus members (Bertout & Genova 2006). Indeed, Dolan 
& Mathieu (2001) found in their kinematic study of A 
Orionis a grouping of candidate members with strong Li 
absorption and u ra d values similar to that of Taurus, al¬ 
though with random spatial distribution. However, Tau¬ 
rus is much closer than A Orionis at a distance of ~ 140 
pc. Moving HD 245185 ~ 3x closer would prescribe 
a stellar luminosity several factors below the main se¬ 
quence for the well-determined A0/A1 spectral type of 
HD 245185 (e.g., Manoj et al. 2006; Folsom et al. 2012) 
or require ~ 2.5 magnitudes of extinction. 


Age estimates for HD 245185 are consistent with the 
~ 5 Myr age of A Orionis, but also within the range of 
ages seen in Taurus (Andrews et al. 2013). Several au¬ 
thors have estimated the age of HD 245185, for example: 
6.9 ±2.5 (Alecian et al. 2013) and 5.5 ±2.0 (Folsom et al. 
2012). Moreover, Franciosini & Sacco (201l) showed that 
HD 245185 sits on the A Orionis cluster sequence in op¬ 
tical and IR color-magnitude diagrams and color-color 
diagrams (see their Figure 2). 

Based on the above information, we cannot definitively 
say whether HD 245185 is a true member of the A Orionis 
cluster, or if it is an interloper from a nearby star-forming 
region such as Taurus. 

4.4. Reconciling HD 245185 with planet formation? 

Disk simulations predict that forming planets clear 
their orbits (e.g., Zhu et al. 2011) producing SEDs with 
rnid-IR dips corresponding to resolved cavities in sub- 
mm/mm images (Andrews et al. 2011); this is obser- 
vationally supported by protoplanetary candidates im¬ 
aged inside transition disk cavities (e.g., Kraus & Ireland 
2012). Dust is thought to be filtered through the cavi¬ 
ties, such that large grains (> sub-mm) are built up at 
the outer cavity edges while smaller particles and gas are 
allowed to flow across (e.g., Rice et al. 2006). 

These signatures of disk clearing from planet formation 
may explain our observations of HD 245185. The mid-IR 
dip in the SED (Figure 5) would originate from planet(s) 
clearing cavities in the disk. The high observed dust 
mass (Section 4.1) would result from sub-mm/mm grains 
containing the bulk of the dust mass being trapped in 
the outer disk, where they are most effectively probed 
by sub-mm surveys. The gas depletion (Section 4.3.2) 
yet ongoing accretion (Section 4.3.3) would be from dust 
filtration allowing gas (and small particles) to continue 
flowing onto the star. 

5. SUMMARY & FUTURE WORK 

We presented results from an 850 -pm survey of the ~ 5 
Myr old A Orionis cluster using JCMT/SCUBA-2. Our 
~ 0? 5-diameter survey region, with a mean rms noise 
level of ~ 2.9 mJy beam -1 , contained 36 (out of 59) 
cluster members with IR excesses indicative of circum¬ 
stellar disks. We found only one significant detection, 
the Herbig Ae star HD 245185, with an inferred dust 
mass of ~ 150 M®. The 3cr upper limit on the aver¬ 
age dust mass of the individually non-detected sources 
was only ~ 3 M®. Thus, despite ~ 15 % of A Orio¬ 
nis cluster members bearing circumstellar disks identi¬ 
fied via their IR excesses, the average disk at ~ 5 Myr 
no longer stores sufficient dust and gas reservoirs to form 
giant planets and perhaps even super Earths. Compar¬ 
ing the dust mass distribution of disks in A Orionis to 
those of other star-forming regions with ages spanning 
~ 1-10 Myr showed that the frequency of disks capable 
of forming systems with multiple giant planets drops by 
roughly an order of magnitude by ~ 5 Myr. 

Our survey showed that HD 245185 is a clear outlier 
in the A Orionis cluster, as its sub-mm emission was far 
above our survey sensitivity and we did not detect any 
other disks in the region. Our follow-up SMA observa¬ 
tions of HD 245185 provided upper limits on optically 
thin CO isotopologue lines luminosities, which we used 
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to infer a gas mass of < 1 Mj up . This implies an ex¬ 
ceptionally low gas-to-dust ratio of < 2, which is > 50x 
lower than the presumed initial ISM value. This over¬ 
all gas depletion sets HD 245185 apart from other Her- 
big Ae disks of similar age, which have been found to 
be gas rich. Indeed, although HD 245185’s strong sub- 
mm continuum emission and large Spitzer IR excesses 
pointed to a massive primordial disk, the mid-IR dip in 
its SED as well as its exceptionally low gas-to-dust ratio 
are clear indicators of a more evolved disk. However, its 
dust mass is still orders of magnitude higher than that of 
debris disks. Its optical emission lines indicative of on¬ 
going accretion, despite significant gas depletion, suggest 
we may be catching HD 245185 during the final stages 
of disk clearing. One possible explanation for our ob¬ 
servations is that the disk around HD 245185 is forming 
gap-opening planets. 

The higher sensitivity of the Atacama Large Mil¬ 
limeter Array (ALMA) could possibly determine the 
residual amounts of gas and dust in these older disks 
and thus place stronger constraints on planet-formation 
timescales. ALMA may also spatially resolve the HD 
245185 disk, which could reveal signatures of ongoing 
planet formation hypothesized in this work. 
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